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Quantum spin liquids are highly entangled ground states of quantum systems with an emergent 
gauge structure, fractionalized spinon excitations, and other unusual properties. While these features 
clearly distinguish quantum spin liquids from conventional, mean-field-like states at zero temperature 
(T), their status at T > is less clear. Strictly speaking, it is known that most quantum spin liquids 
lose their identity at non-zero temperature, being in that case adiabatically transformable into a 
trivial paramagnet. This is the case for the U(l) quantum spin liquid states recently proposed 
to occur in the quantum spin ice pyrochlores. Here we propose, however, that in practical terms, 
the latter quantum spin liquids can be regarded as phases distinct from the high temperature 
paramagnet. Through a combination of gauge mean field theory calculations and physical reasoning, 
we argue that these systems sustain both quantum spin liquid and thermal spin liquid phases, 
dominated by quantum fluctuations and entropy, respectively These phases are separated by a first 
order "thermal confinement" transition such that, for temperatures below the transition, spinons 
and emergent photons are coherently propagating excitations, and above it the dynamics is classical. 
Even for parameters for which the ground state is magnetically ordered and not a quantum spin 
liquid, this strong first order transition occurs, pre-empting conventional Landau-type criticality. 
We argue that this picture explains the anomalously low temperature phase transition observed in 
the quantum spin ice material Yb2Ti2C>7. 



I. INTRODUCTION 

Since the discovery of the quantum Hall effect, it has 
been recognized that phases of matter at zero tempera- 
ture can be distinguished by means other than symmetry, 
e.g. in that case by a quantized Hall conductance, related 
to a topological invariant.^ In recent years, the quantum 
Hall state has been placed into a much broader class of 
phases with topological order of various types, or more 
generally quantum order or long range entanglement P 
Theoretically, many such states may be realized as Quan- 
tum Spin Liquids (QSLs): ground states of frustrated 
quantum magnets with long range entanglement (and 
usually, though not necessarily, the absence of symme- 
try breaking) P The key feature of such QSL states is 
that they support non-trivial excitations which cannot 
be created individually by any local operator, and with 
mutual statistics (and often quantum numbers) differ- 
ent from that of the bare electron and composites made 
from it. In some cases, a quantitative measure of the long 
range nature of the entanglement can be devised through 
a study of the entanglement entrop y, a non-local quantity 
defined in terms of the ground stateP^The entanglement 
entropy is very difficult to measure experimentally, how- 
ever. 

While zero temperature (T) phases can be character- 
ized and distinguished by entanglement properties and by 
their excitations, these criteria fail at nonzero tempera- 
ture, the former for obvious reasons. The excitation crite- 
ria is also invalid at T > 0, where the description of the 
system involves the thermal superposition of all eigen- 
states of the Hamiltonian. In the thermodynamic limit, 
even if the density of elementary excitations (provided 
they can even be defined) is small and finite, their number 



is infinite, and the description of an "excitation above the 
thermal ground state" makes little sense. For a system 
with a QSL ground state, which lacks a symmetry-based 
characterization, one is tempted to conclude that at any 
non-zero temperature, there is no qualitative distinction 
of the physics from that of a trivial paramagnet, and con- 
sequently no phase transition on heating from absolute 
zero to high temperature. 

The true answer is more complex. Some understand- 
ing may be gleaned from work on lattice gauge theory, 
a natural framework for QSL states. It is known that, 
in some but not all cases, the low and high temperature 
phases can be qualitatively distinguished, so that a phase 
transition at non-zero T is inevitable. For example it is 
known that the Coulomb U{1) spin liquid in three spatial 
dimensions and the Z 2 spin liquid in two spatial dimen- 
sions need not go through a transition as the temperature 
is increased, while the Z 2 phase in three spatial dimen- 
sions at infinitesimal temperature is not qualitatively the 
same as the paramagnetic phaseP This difference can be 
understood from a picture of the topological defects of 
the QSL states, whose proliferation coincides with and is 
necessary for the complete destruction of the spin liquid 
state. In the 2d Z 2 and 3d U(l) QSLs, the topological 
defects are point-like objects, Z 2 vortices or "visons" in 
the former case and magnetic/electric monopoles in the 
latter case. The pointlike defects are always created with 
non-zero density at T > 0, making these states smoothly 
connected to a paramagnet. By contrast, in the 3d Z 2 
QSL, the defects are Z 2 vortex lines. In the form of small 
loops, such lines do not disrupt the QSL state. The latter 
is therefore only destroyed when infinitely long loops pro- 
liferate, but, due to their infinite total energy (non-zero 
energy per unit length) , they do so only above a non-zero 
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critical temperature. 

The above considerations give only some indication of 
the presence or absence of a phase transition but do not 
prejudge of the existence of conventional first order tran- 
sitions in cases where one is not required. This is a non- 
universal question, a type often frowned upon by theo- 
rists. However, this non-universality has the virtue that 
the answer may shed light on the microscopic physics of 
the system. In this paper, we address it in the specific 
case of the general model for nearest-neighbor quantum 
pyrochlore antiferromagnets studied in Refs l7l8l which 
contains a U(l) QSL phase at zero temperature, and 
which can be considered a concrete model for "quantum 
spin ice" . The model and the question of a thermally 
driven phase transition is particularly relevant to the case 
of Yb2Ti2C>7, which appears in some experiments to ex- 
hibit a QSL ground state and also clearly shows a sharp 
phase transition in the best quality samples PIUl While on 
the general grounds mentioned above there is no require- 
ment for a T > phase transition, and indeed there is 
none in the simple U(\) lattice gauge theories, we argue 
that the quantum spin ice model model does indeed ex- 
hibit a first order phase transition. The difference is that 
quantum spin ice is described by a U{\) gauge theory 
strongly coupled to "spinons" , fractional spin excitations 
which appear in the theory as matter fields carrying the 
U(l) gauge charge. Thus in quantum spin ice the matter 
matters. 

The first order transition from the QSL to the high 
temperature phase occurs without any change of sym- 
metry, and is thus an analog of the liquid-gas transition. 
Indeed, pushing this analogy, we argue that the transition 
may be regarded as a catastrophic collapse of the QSL 
state occuring at T < T c , which is supported by quan- 
tum coherence, to a thermal spin liquid state for T > T c , 
supported instead by a large residual entropy. In fact, 
the tendency of the thermal spin liquid to supplant the 
quantum coherent phases is so strong that the collapse 
transition persists even in regions where the ground state 
is not a QSL, but an ordered ferromagnet or antiferro- 
magnet (which appear as Higgs phases in the theory). 
We will return to the implications of this finding for the 
putative Higgs transition observed in Yb2Ti2C>7 at the 
end of this paper. 

Our results are summarized by the three-dimensional 
phase diagram shown in Fig. [T] which includes two 
exotic phases, namely a U(l) Coulomb quantum spin 
liquid (QSL) and a U(l) "Coulombic ferromagnet" 
(CFM), wh ose p roperties are now well known at zero 
temperature!^^ Our calculations extend the zero tem- 
perature diagram in the J±/J zz — J z ±/J zz plane to in- 
clude the temperature axis, T/J zz , and as advertised 
above this diagram contains a large-entropy thermal spin 
liquid phase (TSL), in addition to the zero temperature 
like phases. The boundaries in Fig. [I] are calculated us- 
ing the extension of gauge Mean Field Theory (gMFT) , 
described below. While one may be concerned about 
possible artifacts due to this approximation, we provide 



physical arguments that phase boundaries we obtain are 
qualitatively correct. Notably, the transition to the TSL 
is first order and occurs at a temperature strikingly lower 
than the natural energy scales such as the exchange cou- 
plings themselves and the Curie- Weiss temperature. In- 
deed, in the perturbative regime with J± <C J zz and 
J z ± = studied by Hermele et al^ analytic arguments 
imply it occurs at k B T c ~ J±/J zz < J Z z- The gMFT ap- 
proximation actually overestimates T c in this limit, giv- 
ing fcsr| MFT ~ J±/J zz , but does qualitatively capture 
its smallness relative to natural energy scales. 
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FIG. 1: (Color online) Finite temperature gauge mean field 
phase diagram obtained for J±± — and J zz > 0. "QSL" , 
"CFM", "FM", "AFM" and "TSL" denote the (7(1) Quantum 
Spin Liquid, Coulomb Ferromagnet, standard ferromagnet, 
standard antiferromagnet, and Thermal Spin Liquid, respec- 
tively. The lines represent fixed- J±/J zz cuts and are high- 
lighted here to serve as guides to the eye. Details of how this 
figure was obtained are given in Appendix [F] 

We proceed as follows. We first set up T > gMFT, 
mostly extending the analysis introduced in Ref0 and 
present our results after describing the methods used to 
obtain the three-dimensional phase diagram of Figure [T] 
Finally, we discuss our results in the context of the QSL 
candidate Yb2Ti207, whose Hamiltonian is known quan- 
titatively. 

II. GAUGE THEORY 

A. Formulation 

In this section, we recapitulate the spin Hamiltonian 
and its exact slave particle reformulation introduced in 
Refill The Hamiltonian of the system is 

H = J2 [jz.SfS'j - J±(S+SJ + SrS+) 

(ij) 

+ J±± [7«S+S+ + 7 ^SrST] 

+ J^[S?(C«St + £ST)+foj]], (!) 

where the sans serif characters denote components of 
the spins in the local pyrochlore bases, where 7 is a 4 x 4 
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FIG. 2: (Color online) Cuts through the three-dimensional finite-temperature gauge mean field phase diagram obtained for 
J±± = and J zz > 0. "QSL", "CFM", "FM", "AFM" and "TSL" denote the (7(1) Quantum Spin Liquid, Coulomb Ferro- 
magnet, standard ferromagnet, standard antiferromagnet, and Thermal Spin Liquid, respectively. Subfigure (a) (resp. (b), (c)) 
shows cuts for fixed values of T/J zz (resp. J z ±/ J zz , J±/J zz ). Details of how this figure was obtained are given in Appendix [F| 



matrix with only non-zero off-diagonal entries, wh ich are 
complex unimodular numbers, and £y = — T^p^ whose 
explicit expression as well as those of the local bases used 
m Eq. Q are given m Appendix [Aj Details of the gMFT 
formalism were given in Ref |T| and we give here the main 
definitions and results. The gauge "charge" on each dia- 
mond site is 

Qr = Vr S r , r +r) r e„i (2) 

where r/ r = 1 (resp. —1) for a I (resp. II) diamond 
sublattice site, and the are the four nearest-neighbors 
of an rj r (I) diamond sublattice site. The Hilbert space 
is enlarged, and the spins are rewritten 

Sj^r+e^ = S i\r+e (l ^r+e,, j SJ jr+e(j = S r ^ r+e ^ , (3) 

which will allow implementation of mean field theory, 
while preserving the possibility for the description of ex- 
otic phases. Here $ r = e~ lipT is a bosonic spinon field, 
with [cp r , Q r ] — i. For J±± = 0, within this exact refor- 
mulation the Hamiltonian becomes a nearest- and next- 
nearest-neighbor hopping Hamiltonian for the spinons, 
in a background of fluctuating gauge fields (see Eq. (4) 
of Ref. 7). This Hamiltonian is invariant under the U(l) 
gauge transformation 

for any arbitrary real function r i-> Xr- 



B. Gauge mean field theory at T > 

Ref [7| introduced a "gauge mean field theory" (gMFT) 
at T = to decouple the matter fields from the gauge 
fields. Here we extend this analysis to T > 0. Following 
RefHwe now make the Ansatz 

(sp = ssin^e^, (s~) = scos0, (5) 

where the expectation value of an operator U is (U) — 
TyTr [Ue~ f3H ], with /3 = l/(k B T) the inverse tempera- 
ture for Boltzmann's constant ks and Z = Tie~^ H the 
partition function, /i = 0,..,3 and e = (1,1,-1,-1). 
This Ansatz assumes translational invariance (as seen 
in experiment), and is compatible with FM polarization 
along the (global) x axis ((s+) = (s~)). < s = |(s)| < 
1/2 represents the magnitude of the expectation value of 
the spin, which we know in two limits, s(T = 0) = 1/2 
(as in RefJT]) and s(T = +oo) = 0. The variation of the 
magnitude s (and in particular the possibility for it to be 
zero) is the key new ingredient for T > 0. 

To avoid spurious solutions, here we do not solve the 
consistency equations but rather calculate and minimize 
the variational free energy 

F v = F + {H - H ) , (6) 

where now Fo is defined to be Zq = e - ^ , where 
Zq = Tr e~P H ° is the partition function of a fiducial sys- 
tem (we know that the extrema of the variational free 
energy should be the solutions of the mean field consis- 
tency equations) P3 The expectation value of an operator 
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U with respect to the trial Hamiltonian Hq is defined to 
be as usual (U) = ^Tr \U e~^ H °]. As trial Hamilto- 
nian, we choose a decoupled one made of a sum of a free 



spin Hamiltonian and that of a simple nearest-neighbor 
and next-nearest-neighbor hopping Hamiltonian on the 
diamond lattice, i.e. H = + i? s ° with 



(7) 
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{££(*/• ^ $ r+e (1 + h-C.) + E £ *U*r + ^C.) 1 
I r£l |U rell |U J 



and 



EE h . 



(r) 



(8) 
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where J, £l„, t^, h M are real variational parameters, 
which we constrain below. 

We need to determine the free energy F — — i lnZ . 

The free spin part F s ° is trivial, and F% is obtained as 
described in Ref|7] and Appendix [B] (see in particular 
Eqs. (|BT]) and jB5])). Now, trying J — J zz , t'^ = t' 
J±s 2 cos 2 0, and = e^t = 
tional wavefunction, we get 



s 2 sin 29 in our varia- 



F v /N u . c . = 

2{2T[(i + s)m(| + s) + (| 
1 



(9) 
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= 72J^JA-L k ± |M k 
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where u;^ 

IT c o s ( k ' ( e M ~ e v )), M k = J 2 E M e„e ik ^ with 

Ji = 2J±s 2 cos 2 (9, J 2 = 2J 2± s 2 sin26>, e = (1,1,-1,-1). 
We will also be using L k = \ J^^u^ cos ( k ' ( e p _ e ^))' 
M k = J2 tl £ ^ elk • Like in ReffTj A is a Lagrange pa- 
rameter present to enforce the constraint on the spinons 
(rotor operators) <l> J <E> r = 1, in the form ($J$ r ) = 1, i.e. 



1 = Is = 



where 



J 7 ^ = coth 



1 

2K 




A - G 



A - 



(10) 




£ = Ik T |M, 



kl- 
lll) 



Amin + ^1 T 



In the condensed phases, we find A — , 

where A m i n = max k and S = O(l), with S positive and 
independent of T (the difference in the exponent com- 
pared with the zero temperature case explored in Ref. 7 
is addressed in Appendix [Cl . 



Taking the T — > limit of F v and comparing with the 
ground state energy found at zero temperature is rather 
subtle, and carefully described in Appendix |P| We find 
lim T _> i; = -2A(T = 0) + £ l=±1 u/(T = 0), where 
A(T = 0) is determined by the Ia(T = 0) = 1 equation at 
zero temperature, and which is a variational form of the 
zero-temperature ground state energy (see Appendix |d|. 



III. RESULTS 
A. Phase diagram 

We find the phase diagram presented in Figure [l] It 
contains the "continuation" of the four phases present 
at zero temperature: the conventional ferromagnet and 
antiferromagnet, the deconfined U(X) quantum spin liq- 
uid and Coulomb ferromagnet (CFM). This diagram is 
also enriched by the additional thermal spin liquid (TSL) 
phase mentioned in the introduction, which exists at tem- 
peratures T c < T <C J zz ■ 

Within gMFT, we find that the transition to the TSL 
is first order and, at small J±/J zz and J z ±/J Z z, occurs 
when ksT c ~ J 2 /J 2 Z . More precisely 



T c Jz 



3 cos 6 
16 In 2 



4 



sin 2 29 2 



(12) 



where 9 needs to have been well chosen to minimize the 
T = energy (see Appendix [E|). Using the perturbative 
limit of the theoryP the transition temperature would 
be expected to scale as /csTP crt ~ J|./J 2 Z . This means 
that mean field theory overestimates the magnitude of 
transition temperature. 

The zero-temperature properties of the U(l) spin liq- 
uid and CFM were described at length in Refs l7)8|16|17l 
The elementary excitations in these phases are decon- 
fined fractional particles - the spinons and monopoles 
- as well as a gapless photon, which arises thanks to 
fluctuations of the electric field and vector potential.^ 
Hallmarks of those excitations can in principle be seen in 
inelastic neutron scattering. The former two appear as 
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a diffuse signal, and the photon as a sharp, linearly dis- 
persing mode whose amplitude vanishes on approaching 
k = 0. Within the gMFT approach, despite the fluc- 
tuations of the gauge and "electric" fields, these phases 
are nevertheless described by a nonzero (s~), thanks to 
(<&) = 0. In the mean field sense these phases "survive" 
at low but non-zero temperature, i.e. we retain ($) = 
and (s~) 7^ 0, with the latter expectation value reduced 
in magnitude by thermal fluctuations (in fact the reduc- 
tion is very small for all temperatures below T c ). This 
is consistent with the notion that the topologically non- 
trivial spinon and monopole excitations - those generat- 
ing long-range electric and/or magnetic fields - are dilute 
at low temperature, due to their non-zero energy cost 
(gap). A "black body" spectrum of thermally excited 
artificial photons will also be produced by thermal fluc- 
tuations, but as the photons are themselves weakly inter- 
acting (they interact only via anharmonic terms whose 
effects are small at low energy), the thermally excited 
photons do not induce significant scattering. Physically, 
the neutron structure factor should also remain qualita- 
tively similar to its form at zero temperature, modified 
mainly by small thermal rounding. 

Let us turn to the TSL phase. Here the gMFT solution 
is qualitatively changed, and with (s~) = and likewise 
(h M ) = 0. The former implies that, at the mean field level, 
the spinons cannot hop (recall that at least one factor of 
(s*) enters every spinon hopping amplitude), and the 
latter implies that the spins s are freely fluctuating ther- 
mally. Consequently we can view the TSL state as one in 
which spinons are non-propagating, and where there is a 
large true entropy. Physically, this is the best the mean 
field theory can do to emulate the situation in classical 
spin ice, in which the spins are completely free apart from 
the two-in/two-out constraint. This constraint is itself re- 
laxed slowly as the temperature is raised from well below 
J zz to well above it. Hence the TSL state is adiabatically 
connected to the paramagnetic phase, which is described 
by the same order parameter values within gMFT. In- 
deed, if the T c becomes sufficiently large, i.e. comparable 
to J zz , as it will deep in the FM or AFM regimes, then 
the paramagnetic state becomes trivial, and features like 
pinch points need not arise. 



B. Validity of the gMFT treatment 

Within the gMFT solution, the TSL appears via a 
strong first order transition from the exotic states at 
lower temperature. One may be suspicious of this con- 
clusion, since a mean field treatment of related lattice 
gauge theories sometimes gives spurious first order tran- 
sitions. For example, numerical studies of the simplest 
pure compact U(l) gauge theory without matter fields 
show that it undergoes a smooth evolution from T = + 
to high temperature, without any phase transition. Nev- 
ertheless, a mean field treatment predicts a first order 
transition in that case as well.^ It is therefore natural 



to ask: What is wrong with the mean field treatment in 
this case? Does the same problem affect our calculations 
for the pyrochlore problem? 

In the pure compact U(l) gauge theory, the Hilbert 
space consists of the space of electric field configurations 
on the lattice which strictly satisfy the charge neutrality 
constraint div E = on all sites. The Hamiltonian in this 
case has the usual form, 

h u(1) = uJ2 e ^-^J2 cos ^ 7x A )p> ( 13 ) 

(v) p 

where the second sum is over spatial plaquettes p, and 
gives an explicit microscopic stiffness K penalizing con- 
figurations with non-zero magnetic flux B p = (V x A) p . 
As usual, Aij and Eij are canonically conjugate variables. 
The mean field treatment consists of decoupling this stiff- 
ness term by defining a self-consistent value of the bond 
"order parameter" (e ij4iJ ), calculated using a Hamilto- 
nian for decoupled bonds. At the mean field level, there 
are two phases: one where this expectation value (e %Aii ) 
is non-zero, defining a putative "Coulomb phase" , and 
another where the expectation value vanishes, signaling 
a "confined" phase. This transition occurs abruptly, i.e. 
is first order, because the mean field treatment essentially 
neglects correlations amongst the fluctuations of the elec- 
tric fields. In reality, if one considers the situation with 
large K and low temperature, the electric field fluctu- 
ations are highly correlated. They arise from magnetic 
monopole excitations, i.e. textures in B p , whose mo- 
tion leads to long-range electric field fluctuations through 
Faraday's law. The complicated spatial structure of mag- 
netic monopoles is completely missed in the mean field 
treatment. Instead, in mean field theory, the electric 
field fluctuations occur locally (and do not even obey the 
diverge-free constraint since the latter is only satisfied on 
average in mean field theory). In fact, because monopoles 
have a finite energy, at any T > they appear in a 
non-zero (albeit exponentially small) concentration, and 
their motion immediately leads to the destruction of the 
Coulomb phase, in a strict sense. Specifically, they screen 
the interactions between inserted test magnetic charges, 
and lead to exponential decay of all correlations. Never- 
theless, at very low temperature, because the monopoles 
are very dilute, the short distance stiffness given by K 
is largely unaffected. With increasing temperature, the 
density of monopoles increases, gradually reducing the 
effective stiffness on the scale of the correlation length 
(distance between monopoles). We can conclude that 
mean field theory fails in this case because it misses the 
true mechanism of destruction of the Coulomb phase at 
T > 0, which is magnetic monopole proliferation. 

Now let us turn to the real pyrochlore problem. Here 
the gauge theory is distinguished from the pure t/(l) 
problem described above by the presence of matter 
(spinon) fields $ r (and by the constraint on the magni- 
tude of the electric fields, but this is of secondary impor- 
tance). Most importantly, there is in fact no microscopic 
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stiffness for the magnetic flux. Recall that s ± ~ e ±lA rep- 
resents the gauge magnetic vector potential in this for- 
malism, and these fields appear only in the kinetic terms 
of the spinon variables. Hence in this problem the mag- 
netic stiffness arises only dynamically: spinons propagate 
coherently most efficiently through a background of zero 
magnetic gauge flux, and thereby have lowest kinetic en- 
ergy in that situation. It is the lowering of spinon kinetic 
energy that is responsible for the magnetic stiffness, and 
hence stabilization of Coulombic (QSL and CFM) phases 
at zero temperature. At T > 0, there are now two sources 
of gauge fluctuations, in contrast to the situation in the 
pure gauge theory. Gapped magnetic monopole excita- 
tions still exist (in these phases), but in addition we may 
have thermal activation of excited spinon states. The for- 
mer process is similar to that in the pure gauge theory, 
and is missed by gMFT. The latter process is captured 
by gMFT, and acts to reduce the microscopic magnetic 
stiffness even on short length scales. As this stiffness is 
reduced, electric field fluctuations grow in response, fur- 
ther decreasing the stiffness leading to a rapid explosion 
of gauge fluctuations and rapid reduction of the spinon 
bandwidth. Once it reaches the thermal energy fc^T, 
the Coulomb phase collapses entirely as there is no mi- 
croscopic stiffness to support it. 

In reality, both this process and the one due to ther- 
mally excited monopoles should be responsible for de- 
struction of the Coulomb phase. Since gMFT captures 
one but not both of these mechanisms, we expect it to 
be a better approximation here than in the pure gauge 
theory, but still susceptible to possible O(l) errors. Our 
expectation is that, like most mean field theories, the ne- 
glect of correlated fluctuations will lead to a reduction of 
the true critical temperature in comparison to the mean 
field result, but likely not suppress the transition entirely 
Some further evidence for this conclusion comes from ex- 
amining the phase diagram more broadly. 

First, let us consider the role of symmetry. The above 
discussion, and comparison to the situation in the pure 
gauge theory, applies best when the ground state is in the 
QSL phase, which breaks no symmetries and at T > 
can be adiabatically connected to the paramagnetic state. 
In the other phases, which break symmetries, a phase 
transition is required at T > 0. In crude but physi- 
cal terms, we can imagine two possible scenarios. On 
one hand, confinement may occur at a lower temperature 
than the restoration of symmetry, so that a conventional 
Landau picture describes the symmetry breaking transi- 
tion. On the other hand, confinement may occur simulta- 
neously with symmetry restoration, which is what occurs 
in gMFT. In the former case, we should expect that the 
phase transition should be described approximately by 
the usual Curie- Weiss mean field theory (CWMFT) in 
terms of self-consistent exchange fields and spin expec- 
tation values,^ since once the transition temperature is 
reached, confinement is occuring on short length scales, 
and the microscopic spin variables are good order pa- 
rameters. It is interesting to compare the CWMFT tem- 



perature to that predicted by gMFT. We find that the 
T c from CWMFT is systematically significantly larger 
than that found in gMFT. For example, for the parame- 
ters corresponding to Yb2Ti207, the critical temperature 
in CWMFT is T C CWMFT = 3.2 K,E while in gMFT it is 
T gMFT = o 56 K The much smaller value f the cr itical 

temperature in gMFT is strong evidence that confine- 
ment physics plays a role in the transition. Note also 
that the gMFT value is much closer to the observed one 
in Yb2Ti207, supporting this notion for experiment as 
well. 

Second, we may consider the role of Higgs condensa- 
tion. When J± and/or J z ± are not too small, the ground 
state is not a deconfined but a Higgs phase, with a spinon 
condensate. Consequently, the magnetic stiffness is en- 
hanced beyond the usual dielectric form to a Meissner 
one, such that magnetic gauge flux is actually expelled 
from the system. In particular, in these Higgs phase (the 
AFM and FM in Fig. [I), due to this Meissner effect, 
the energy of a magnetic monopole is no longer finite 
but actually infinite. Consequently, magnetic monopoles 
cannot be thermally activated in these phases. There- 
fore only the mechanism of spinon fluctuations (weak- 
ening of the Higgs condensate), which is captured by 
gMFT, is present in these regimes, and we expect the 
accuracy of gMFT to increase. The fact that a single 
T > confinement-like transition appears here and is 
smoothly connected to the one appearing for small J±, 
J z ± couplings suggests that gMFT is qualitatively cor- 
rect throughout the phase space. 

Based on the above reasoning, we conclude that the 
T > transitions for quantum spin ice in the corner of 
phase space studied here are qualitatively correctly de- 
scribed by the gMFT treatment, and should be thought 
of as confinement or quantum-to-classical transitions. As 
a consequence of the latter, the experimental signatures 
(such as "pinch points" - see Sec. IV A) of classical spin 



ice are expected in the TSL regime above the critical 
point. In particular, the famous "pinch points" should 
appear in this phase. 



IV. DISCUSSION 

A. Connection with real materials: the case of 
Yb 2 Ti 2 7 

The Hamiltonian parameters J zz , J±, J z ± and J±± 
(see Eq. 0) of Yb 2 Ti 2 7 were extracted in Refill by 
fitting linear spin wave theory to high-resolution in- 
elastic neutron scattering in high field. The accuracy 
of the values J zz = 0.17, J± = 0.05, J z± = -0.14 
and J±± = 0.05 meV was subsequently confirmed in 
Refsfiniini through comparison of high-temperature spe- 
cific heat and entropy data (for v arious exchange param- 
eters reported in the literature)! 8 * 11 * 15 ! Despite the evi- 
dent complete quantitative knowledge of its Hamiltonian, 
the nature of the low-temperature phase of Yb2Ti207 



7 



in zero field is still open to debate. Several studies 
find no sign of order down to the lowes t accessible tem- 
peratures (30 mK in Ref HU) EEMmiSl anc l diffuse neu- 
tron scattering at T — 30 mK and H = compati- 
ble with a two-spinon continuum!^ Two other neutron 
scattering studies have repo rted t he presence of an or- 
dered ferromagnetic moment |H12U Specific heat measure- 
ments reveal strong sample dependence, which has re- 
cently been associated to Yb substitution ( "stuffing" ) on 
the Ti sitep^l so it is possible that such disorder modifies 
the zero field ground state in some samples. However, 
even this is not clear. Not knowing for sure what the 
low-temperature phase is, the nature of the transition 
observed at T ~ 200 mK! 9 ' 11 ' 12 ' also remains equivocal. 

Recently, experimental and theoretica l 10 * 19 ! works 
have also addressed the nature of the phase transition. 
The authors of Ref |TT| argue that their experiments im- 
ply an ordered ferromagnetic phase, and provide evidence 
for the first order nature of the thermal transition to this 
phase. In ReffTOl a theoretical model with third neigh- 
bor exchange (which can be considered a perturbative 
approximation to the full H)^ and consequently a ferro- 
magnetic ground state is studied by Monte Carlo sim- 
ulations, finding a first-order transition. In our gMFT 
calculations, as discussed above, all the transitions with 
increasing T are first order as well. 

Ref(TT] suggests that the thermal transition in 
Yb2Ti2C>7 may be regarded as a "Higgs transition". We 
would like to discuss this interpretation, in light of the 
one we have offered above. The term Higgs transition has 
an accepted meaning in quantum field theory, where it 
refers to a transition which may be described as the con- 
densation of a bosonic field carrying a non-zero gauge 
charge, and coupled to a dynamical gauge field. In our 
formulation, such a bosonic field is the spinon, which 
carries the electric gauge charge, and phases with non- 
zero spinon condensates are indeed Higgs phases, and 
correspondingly have magnetic order. The zero temper- 
ature quantum phase transitions from the CFM and QSL 
phases into FM and AFM phases are indeed (quantum) 
Higgs transitions in this sense. In a strict sense, the situ- 
ation at T > prohibits any true Higgs transitions, since 
the Coulomb phase itself is not sharply defined at T > 0, 
i.e. there are no critical gauge fields anywhere in the 
phase diagram at T > since the U(l) gauge fields are 
compact. Even if we look for a non-strict interpretation, 
since the CFM and QSL phases are deconfined, tran- 
sitions from them to the confined paramagnetic phase 
with increasing temperature can clearly not be regarded 
as Higgs transitions, as none of these phases have Higgs 
condensates. A non-strict view of the thermal transition 
from the FM or AFM phases to the paramagnetic one as 
a Higgs transition might be possible. However, since the 
FM and AFM phases are already confined states even at 
T = 0, it seems unreasonable to consider spinon conden- 
sation as the mechanism for this transition. 

From our point of view, the essence of this transition 
is not Higgs condensation but confinement, as discussed 



extensively in the previous section. According to this 
picture, supported by gMFT, the phase above the tran- 
sition should be regarded as a classical thermal spin liq- 
uid (TSL) , similar to low temperature regime of classical 
spin ice, and the transition may be regarded as describ- 
ing the release of entropy associated with an abrupt loss 
of quantum coherent spin dynamics. We would therefore 
expect significant differences in the inelastic spin correla- 
tions below and above the transition. Specifically, above 
the transition diffuse scattering with visible pinch-point- 
like features would be expected, provided T c -C J zz , so 
that the classical ice rules are not strongly violated. Be- 
low the transition, quantum coherence is known to wash 
out the point points, even in the QSL regime. Indeed, fea- 
tures reminiscent of pinch points were observed above T c 
in Ref ITTI in Yb2Ti2C>7. Another piece of support for this 
interpretation is the rough agreement of the critical tem- 
perature predicted by gMFT for this confinement tran- 
sition, 7 1 s MFT « 560 mK, with the experimental value 
ycxp py 255 mKP^The agreement is much better than in 
the usual Curie- Weiss MFT for which T C CWMFT = 3.2 K, 8 
as noted above. 

A more tricky issue is the presence or absence of mag- 
netic order below T c , which is controversial experimen- 
tally, and is likely to be related to the Yb stuffing men- 
tioned earlier. It is not unreasonable to expect the extra 
Yb spins to affect the state and dynamics of the sys- 
tem. Yet, if so, it is surprising that the microscopic 
model which neglects them appears to provide an ex- 
cellent quantitative description of a comprehensive set of 
data in applied magnetic fields and/or higher tempera- 
tures. A possible interpretation is that the "interstitial" 
spins become polarized in modest applied fields, and fluc- 
tuate paramagnetically at higher temperatures, in either 
case obviating their effects on the Yb sublattice spins. At 
low temperature and low fields, however, the interstitial 
spins may be free to couple to the main sublattice, and 
act to effectively modify the exchange parameters, mov- 
ing the material around in the general phase diagram, 
perhaps even into the QSL regime. From Fig. [I] one sees 
that although the ground state is quite sensitive to val- 
ues of the exchange parameters, the thermal transition 
remains confinement-like regardless of these values, con- 
sistent with observations. Obviously this interpretation 
is highly speculative. Disorder might have many other 
unanticipated effects, and understanding the presence or 
absence of magnetic order and the sample dependence of 
various experiments requires significant more experimen- 
tal and theoretical study. 



B. Conclusions 

We have studied the development of a quantum spin 
liquid ground state and its neighboring phases on cool- 
ing from high temperature. We argued that in the case of 
quantum spin ice, quantum coherence of spins onsets in 
an abrupt first order transition. When the ground state 
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is a deconfined spin liquid, this transition may be viewed 
as describing the confinement of fractional spinon excita- 
tions. Above the transition, the system has substantial 
entropy and behaves as a thermal spin liquid, with many 
of the characteristics of classical spin ice. 

In quantum pyrochlore magnets, more experimental 
studies on the spin correlations at intermediate temper- 
atures would be welcome, not only from neutron scatter- 
ing, but also other probes of spin dynamics such as muon 
spin resonance. A thorough study of the development of 
the phase transition in Yb2Ti 2 0y with applied magnetic 
fields would also provide considerable fuel for future the- 
oretical work. The techniques of the present paper, along 
with other approaches, can certainly address such prob- 
lems. 

This raises the more general question of the existence of 
first order quantum to classical transitions for other mod- 
els and materials with spin liquid ground states. Most 
of the techniques used to study quantum spin liquids 
address either ground state properties (e.g. Gutzwiller 
variational wavefunctions, Density Matrix Renormaliza- 
tion Group) or are limited to relatively high temper- 
atures (e.g. high temperature series expansion, quan- 
tum Monte Carlo), leaving the intermediate temperature 
regime where such a transition might occur relatively un- 
studied theoretically. It may consequently be interesting 
to develop theoretical methods for such temperatures in 
the future for frustrated quantum spin models. 
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is a third root of unity. 



(A2) 



Appendix B: Details of the calculations 



Here we describe the calculations leading to Eq. ([9]) in 
great detail, and proceed making simplifying assumption 
as we go. 

The decoupled Hamiltonians H% and H® such that the 
trial Hamiltonian is H° = + H® are given in Eqs. ^ 
and ([8]). From them, we need to determine the free energy 
Fo = —jjlnZo = Fq + Fq. The free spin part is trivial, 

it is the free energy of free spins in a field h: 
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-T 



2 cosh 



(Bl) 



ln2--ln(l-4|(s- r , r+e J| 2 ) 



Appendix A: Notations 

The local cubic bases in which the Hamiltonian Eq. ([!]) 
is expressed are the following (Hj, b,-,ej) bases 



e = (1,1,1)/V3 
ei = (1,-1,-1)/V3 
e 2 = (-1,1,-1) A/3 
e 3 = (-l,-l,l)A/3, 



fio = (-2,l,l)/V6 
ai = (-2,-l,-l)/V6 
&2 = (2,1,-1)/V6 
a 3 = (2, -1, l)A/6 

(Al) 



To compute F*, we rewrite the first part of H^, as 7 



I, II r£l,II ^ ' 



re 



where the second part of the right-hand-side is intro- 
duced to implement (a relaxed version of) the constraint 
$t$ r = 1 ( S ee Refill). Therefore, A 6 K + serves as a 
Lagrange multiplier. We get 
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|^ntn r + A ($t$ r -1)U E E ^ Sj+e^r+e, + E E C ^-e^r-e„ | (B3) 
- I E E (*/* ^ $ r+e tl + h.C.) + E E (> ^r-e^r + h.C.) 1 , 



Then, using the following Fourier transformation conven- 
tions 



x k 



— £■ 



,t 



rrl>l< 



x k e 



ik-r,- t 



E4e- ik ' r S (B4) 



where fL 



is the bosonic Matsubara frequency, we 



find the Green's function of Hq, 



[Go 1 ] = 



ij^l + A - Li, 



-Ml 



1 , , 2 



-M k 

A-i k , 



, (BIO) 



we arrive at 



2AiV„ 



k i=± 



i=± k 



(B5) 

where w k involve i and t' introduced in Eq. ([7]), and J, 
and where the last term comes from the constant term in 
Eq. (B2 1. Note that A and are quantities defined from 
the fiducial system with Hamiltonian Hq (in Ref|7] their 
"equivalents" were defined from the gMFT Hamiltonian) . 



where 
M k 



2 EV* k ' e "> (BU) 



setting t* = t Vfl . Here L and M are defined in terms of t 
Let us now calculate (H) and (H Q ) . Going to the and t', which are at this stage arbitrary. In the main text 



action formalism, with the conventions 

"+P/2 



xn„ 



1 

2^ 
1 

2^ 



dr x T e 



-13/2 
+13/2 

-13/2 

X n 6 



we give formulae for these quantities which are equal to 
the above ones when the specific values for t and t' have 
been taken. For simplicity we use the same symbols for 

(B6) 

both expressions. From this we get the spinon dispersion 
relations 



(B7) 
(B8) 



u>± = V2J^\-L k ±\M k \, (B12) 
Note that, like A and oj, L and M are defined here from 



(B9) the fiducial Hamiltonian. So, inverting Eq. (BIO I: 



G (k, 00 



in 2 



X-L, 



. / 



Ml 



2.7 n 



+ A - L k , 



|M k | 



-1 



1 



-1 



i i 



1 



(B13) 



(BU) 



where 



Carrying out the sum over all integers n, we get 



n 2 

2.7 
ik T |M k 



A - Li 



- |M k | 2 (B15) 
(B16) 
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Go(k,T = 0) = ^Go(k,f2 n ) = i 



\\M k \ 



M k 
|M k | 



^,7 



(B17) 



where 



where Go(k) = Go(k, t = 0), so 



J" k = coth 



x-et 



(B18) 



As usual, we have 



<4kVk>o - G?i,i(k). (B19) 



r£l,II I r£l (-i^^f-i r£ll fi,is^f-i 

~ { E E (*M (*J $ r+e> + h.C.) + E E (> (*L, M *r>0 + h.C.) } . 



(B20) 



Now, recall: 



(H)o = ^ E (n^r)" "4 EE^ <*U,*r+e„>0 +E E ^ (^-, M *r-,„>0 l (B21) 
r£l,II I r^I \x,v^[i r£ll \i,v^\x J 

j *±|e E (7;,s^<*i*r+e>+h.c.) + E E (7p;^($u$ r >„+h.c.) , 

r£l il,v^\l r£ll \l,v^fc\i J 



where we have defined x = (x)o, and with 



1 /Jv- J? 



ik-(e M -e,/) 



(B22) 



and 



(^ r+e Jo = EKk $ n,k)oe- 



Jy^MZ 



2N U . C . V 2 ^ |M k , 



! k 



(B23) 
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Finally, 

(H - H Q ) 



(H) -{H*) -{H S } 

J zz J 



E < n ^>° 



(B24) 
(B25) 



rGlJI 



l E E ( J ±VC - O (<+e t . < fr+ej0 + E E ( J ±^ - V) (Ce.VeJo i 



EE 



h.c. 



(B26) 



+ EEM r )-(Sr, 



EE 

r£ll H,v 
r+e„)o- 



Now, simply assuming | (s r ,r+e M ) o I = s to be independent 
of (j,, we have already: 

( ff s°)0 = - E E H ^( r ) • &.r+e> (B27) 



4JV„.c.s. l + 2s 
In 



/9 1 - 2s ' 
and if we use the Ansatz Eq. ([5| , 

= s En sin 9 and s~ = s cos 9, 



(B28) 



then 



E ^±7^s^s* = 2s 2 J z± e l/ cos6'sin6l = s 2 J z± e y sin 20, 



and J±s M s* = J±s 2 cos 2 ' 
We now take 



t' nv = t' and t v = e v t, 



so that Eq. (B24| becomes 



(B29) 



(B30) 



(H - H ), 



J zz J 



E ( n i n -)o 



rGLII 

f -*'J \E E (*UWo+E E <*U*^>o 

r£l fi,u^fi r Gil fi^^fi 



■h 



- 1 



{ E E $ r +e jo + h.c] + E E [m $ Ua>o + h. 



EE^W • (Sr,r+e f ,)o, 



(B31) 



where 



Appendix C: Explicit expression of the 73 = 1 
constraint in the condensed and deconflned phases 



Ji = 2 J ± s z cos 2 9 and J 2 = 2 J z ±s z sin 20. 

(B32) 

Finally, setting J = J Z2 , i' = t' = ^ = J±s 2 cos 2 9 and 



tu^e^t = e u 4y = J z ±s 2 sin 29, and using Eqs. (B22|, 
( |B23[ ) and ( |B28| ) we recover Eq. j^. 



The constraint on the spinons (rotor operators) 
$t$ r = 1 is enforced in the form ($J$ r ) = 1, i.e. 



l = /a 



1 J 



2AT U . C . 



E 



Ft 



A - £,1 



X-£Z 



(CI) 
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where 



with 



coth 



L k T |Af k | 



(C2) 



While in the deconfined phases the sum in Eqs. ( 10 ) 



and (CI) can be turned simply into an integral (since 
A — £^ > for all k), in the condensed phases, the spinon 
dispersion relation hits zero at a wavevector ko, i.e. such 
that £^ o = A m ; n = maxk^k, and one should allow for a 

subextensive part in A. We write A = A m j n + ST/N U . C ,, 
where 8 = 0(1), with 8 positive and independent of 
T. This leads to (see also the Supplemental Material 
ofRefE]): 



ia — +1 



4 



(C3) 



where 1% is the right-hand-side of Eq. (CI) turned into 
an integral and evaluated at A 



A, 




A ■ -r 



and /™ ln is the part of 

does not vanish when N u c 
find, in terms of 8, 



oo. Defining p -- 



i for(?,s,J z± ^0 
1 otherwise 



(C4) 



which 



we 



(C5) 



or more generally, considering the rotor constraint, p — 

Note: (i) the difference in the exponent compared with 
the zero temperature case explored in Refi7:, which comes 
from the coefficient modification involved with J^i fa) 
we used a hat on 8 because its definition differed from its 
"equivalent" at T = 0. 

A final remark is in order: in 
EkE, = ±i [ 2rln hzA of Eq. <9), the k term 

L 1 — e 1 k J 1 — 1 

is of order hxN u , c , (for N u , c , large at fixed T) in the 
condensed phases. Since lniV u . c . = o(N u , c ,), this con- 
tribution is actually negligible at large iV„. c . compared 
with the main contribution to the free energy, which is 
extensive. 



-2I 2 (6, A) cos 2 9J ± - 4/1(6*, A) sin 26 J Z ^D2) 



1 



( !■■ i 



where 



(D3) 



A will have been determined by 
solving J3 = 1, and 6 either by solving the consistency 
equations (and choosing the lowest energy solution), or 
by minimizing the Eq. (Dl ) form of the energy. 



The energy can also be derived through another pro- 
cedure, which we call "by decomposition." The energy 
found in such a way is not variational (i.e. the ground 
state energy cannot be found by minimizing it over 9); 
the ground state energy is found by plugging in the pa- 
rameter values found by solving the consistency equa- 
tions. For those values the decomposition and the vari- 
ational forms of the energy yield equal values. The en- 
ergy is found by analyzing what the decoupled problem is 
equivalent to. The mean-field Hamiltonian, found using 
the usual decomposition (see Eq. (6) of Refill , is 



so that 



We find 



H 



MF 



ttMF 
spinon 



ttMF 
spin 



E. 



MF 



(D4) 



77.MF _ N t Mr , Mr _ Mr \ lj^r\ 

& — ji ii.c. V spinon ' fc spin fc constj ■ \ lJO ) 



MF _ J 
e spinon — ^ j \~k 



2A 



(D6) 



Oi„ : -\ Kv + h MF (D7) 
4 V onst = ~4/ 2 cos 2 6J ± - 8/1 sin 26 J z± , (D8) 



where h, I2, and h^F were defined in RefJT] It turns out 
that we find, empirically (i.e. numerically), 



ye 



MF 



= E, 



GS 7 



(D9) 



provided A is chosen such that J3 = 1 is satisfied. 

Note that, to avoid any confusion, we refrained from 



calling e^f non 



v resp. 



MF ' 
spiny 



(resp. 



as they do 



not come from H^, and H® of the present paper. 



Appendix D: Comparison with zero temperature 

Here we outline the procedure described in Refill i- e - 
in the case of zero temperature, and show that this case 
is recovered when we take T — > in the present work. 



1. Energy at T = as derived in Savary and Balents 



In Ref(7J we reported 

(H)gMFT,T=0 = Eqs 



(Dl) 



2. T — > limit of the variational free energy 

We now go back to the variational free energy at T > 0, 
Eq. (pj), and take its T -> limit. Since lim T -H.o s = 1/2, 
the first part of the first term of Eq. ([9| vanishes in the 
T -> limit. Therefore, 



lim F v 

T-s-0 



-2N U c lim A 



(D10) 



k i=±l 



lim 

T->0 



2Thv 
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Since A is determined through the 1% = 1 constraint, let 
us address the becoming of the latter. It is 



1 = 




(Dll) 



with T defined in Eq. ( 11 ), 



= coth 



(D12) 



In the deconfined phases, since A - ^ > for all k, 
J 7 ^ — > 1 trivially, and it is obvious the I3 = 1 constraint 
reduces to that found at zero temperature in Ref. 7 thanks 
to the choice (described in Section pi B and Appendix |b|) 
of parameters t, if and J of the fiducial Hamiltonian. It 
follows immediately that lim T _, A dccon = A docon (T = 0). 
In the condensed phases, as described in RefJT] and Ap- 
pendix O the sum is better split into a ko term (which 
we call p) for which £^ o = A m ; n = maxk £± > ano - remain- 
ing terms 1$. Defining A = A m i n + ST/N U , C , in the con- 
densed phases and taking the N u _ c , — ¥ 00 limit before 
the T — > limit (since physically we are interested in 
low but non-zero temperature but thermodynamically 
large systems), we find lim-r-^o A m i n = A m i n (T = 0) 
limT-s-o I3 = T 3 (T = 0), which, from p = 1 — I' 3 , implies 
liniT_>o p = p(T = 0), and as a consequence, 



lim A 

T->0 



A(T = 0), and lim w 



■XT 



Finally, the very last term limr->o 



2Tln 



0). 

(D13) 
of 



Eq. (D10) goes trivially to zero in the deconfined phases. 



In the condensed phases, as noted in Appendix [Cj the 
ko term goes as lniV u . c . for N u c large (at fixed T), but 
the zero-temperature limit takes this term to zero so that 



2Tln- 



lim T ^o 
well. 

Finally, we arrive at 



in the condensed phases as 



lim F v 



= -2N U . C X(T = 0)+Y.J2 ^ T = °) 



k i=±l 



(D14) 



Eq. (11) (of Ref; 7), and selected those with lowest en- 
ergy. 

In Ref 13 the consistency equation Eq. (11) as well as 
the variational form of the energy Eqs. (12) and (13) 
(Eqs. (D2) and (D3| here) involve the sums I\ and I2 
which are subject to greater numerical errors than the 
variational form (of the type of -E^^ on ) presented here. 
This led to a small mistake in the position of the AFM- 
FM phase boundary in RefJTl We believe the slightly 
modified diagram presented (Figure [3]) here is correct. 
An Erratum with this correction is being simultaneously 
submitted to Physical Review Letters. 




FIG. 3: (Color online) Zero-temperature gauge mean field 
phase diagram obtained for J±± — and J zz > 0. "QSL" , 
"CFM" , "FM", "AFM" denote the 17(1) Quantum Spin Liq- 
uid, Coulomb Ferromagnet, standard ferromagnet, and stan- 
dard antiferromagnet, respectively. 



Appendix E: Analytical phase transitions in the 
small parameter regime J±,J Z ± <SC J Z z 

Here we look at simple limits and find approximate 
phase transitions analytically, taking advantage of the 
fact that the phase transitions to the TSL occur at low 
temperature. 



Using the empirical evidence Eq. ( D9 1 , this proves that 



the variational form of the ground state energy is re- 
covered when we take the zero temperature limit of our 
variational free energy. Note that Eq. (D14) could also 
be seen as a convoluted proof of Eq. (|D9|)! 



3. Discrepancy with the T — phase diagram 
computed in Savary and Balents 

To obtain the phase diagram of Figure 1 of Refill we 
did not minimize the variational ground state energy 
at T = 0, but rather solved the consistency equations 



For s = 0, 
F„(s = 0)/JV u ., 



and 



1 = 



1. s = 



-AT I In 2 - In 



(El) 



1 — e T 



coth 



1 J ZZ X 
TV 2 



(E2) 
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For small enough temperature, we find A = J zz /2, which 
leads to 



and at 



i^(s = 0)/JV u . c . a J„-4Tln2. 



1/2 



(E3) 



When s = 1/2 (which encompasses but is not a priori 
restricted to the T = case), 



F„(s = l/2)/JV„. c 



(E4) 



F„(T = 0)/iV u . c . - 2T V / In ^ 



so, for small enough temperature, the variational free 
energy at T ^ is almost constant and equal to that at 



T = 0. Eq. (E4) is valid both in the condensed and un- 
condcnsed phases. The absence of a distinct correction 
to this form in the presence of a condensate may be inter- 
preted physically as the fact that the condensate carries 



zero entropy. Since the final term in Eq. ( E4 1 is a mea- 



sure of the entropy, it is not corrected by a condensate. 



T T — 

-L C° ZZ 



In 2 16J Z C ± 2 - 3J| 2 



(E8) 



V3 



for J z ± > — i=J± where 8^0. 
2v 2 

Appendix F: Calculation of the phase diagram and 
representation of 3D surfaces and cuts 



The phase diagram, Figure [T] and the cuts, Figure [2] 
was obtained by sampling points separated by 0.1 incre- 
ments in the J±/J zz direction, 0.0125 to up to 0.05 in 
the J z ±/J zz direction and 0.0025 in the T/J zz direction 
in regions surrounding a phase transition. 

Minimization was obtained by comparing values of F v 
for values s = 0.0001, 0.005, 0.01, 0.02, 0.05, 0.1, 0.2, 
0.3, 0.4, 0.45, 0.46, 0.47, 0.48, 0.49, 0.495, 0.496, 0.497, 
0.498, 0.499, 0.49999, and 6 = rwr/32, n = 0, .., 7 for each 
of the sampled points. Which phase each point belonged 
to was determined according to Table |TJ 



Transition to the s = state (TSL) (if it indeed 
occurs at small T and small J±) 



We apply the results Eqs. (E3| and (E4| to hnd ana- 
lytical forms of the transitions. 



For J z ± = 



16J 22 ln2 



T c J zz \n2 



(E5) 



TABLE I: Criteria for determining the ground state phase. 
Note that, in practice, the criterion for 9 or s zero or nonzero 
is 6 or s smaller or greater than 10 -6 . 



p 


e 


s 


phase 










QSL 







7^0 


CFM 








FM 







7^0 


AFM 











TSL 



b. For J z ± / 

A set £ of phase transition points was subsequently ob- 
tained by taking the midpoints (along well-chosen lines) 
between sampled points not belonging to the same phase. 

T C J ZZ — J± 2 H Jz± 2 i (E6) The phase transition surfaces were then obtained by 

''' - -1 In - triangulating the projections of the points of £ onto ap- 



where 9 needs to have been well chosen to minimize the propriate planes. 

~~ oj> ' ' The triangle edges were subsequently parametrized so 



that any surface cuts could be obtained. 

Lines were 
(E7) polynomials 



3J£ 2 

T c J Z z = — ~z~ — : — - for J z + < — = J+ where 9 = 0, T . . rii 

16J 2Z In 2 2y2 Lines were smoothcncd by fitting cuts to fourth-order 
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24 Throughout the paper, we use the usual vocabulary of 
quantum field theory, which is also that used in Ref ll7l 
This means in particular that our "electric" fields are the 
local z-components of the spins (called "magnetic" fields 
in the classical spin ice literature), and that the (classical 
analogues of the) spinons denote what are referred to as 
"magnetic monopoles" in the classical spin ice literature. 



